Rayleigh's collapsing time of a spherical cavity: the A-factor 

J. A. S. LimJ] 

Universidade de Sao Paulo, Instituto de Astronomia, Geofisica e Ciencias Atmosfericas 
Rua do Matao, 1226, CEP 05508-900, Sao Paulo, SP, Brazil 

?v n 

^ ! F. E. M. da Silveir^ll 

[T j | Universidade de Sao Paulo, Instituto de Fisica 

Caixa Postal 66318, CEP 05315-970, Sao Paulo, SP, Brazil 
i ,! (Dated: February 2, 2008) 

Abstract 

New corrections to the equation of motion and total collapsing time of an empty spherical 

« : 

cavity immersed in an infinite incompressible medium are proposed on the assumption of a non- 

O ■ 

uniform density. The dimensionless number quantifying the corrections with respect to the stan- 

>v 

dard Rayleigh results (coined the A-factor) is fully independent of other possible contributions 
like surface tension and viscous terms. The A-factor effect advocated here can be seen as a direct 
consequence of a mass-shell non-trivial solution to the continuity equation. The consistency of the 
corrections with respect to the Bernoulli theorem and some physical consequences in the framework 

t-H \ 

CS| , of the Rayleigh-Plesset equation are also discussed. 
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I. INTRODUCTION 



Ninety years ago, Lord Rayleigh published his seminal article on the problem of deter- 
mining the collapsing time of an empty spherical cavity, starting from an arbitrary radius a, 
suddenly formed in the bulk of an infinite medium [l| . Later on, that paper became the pre- 
decessor of a very active area of research, collectively referred to as bubble dynamics, and, as 
such, it can be regarded as a kind of paradigm to a vast class of spherically symmetric one- 
dimensional non-steady flows. To quote some examples, Rayleigh's work is closely related 
to a large number of applications in the development of engineering devices I2J, sonolumi- 
nescence js], and, more generally, with the important problem of cavitation jj Nowadays, 
that classical approach is the starting point for a large number of investigations, based on 
the medium properties therein neglected, like surface tension, fluid viscosity, heat transfer, 
acoustic cavitation, bubble relaxation, diffusive terms, shock-wave induced collapse, chaos, 
and many others [5-9] (see also 1(3, 11] for reviews on such subjects). 

Rayleigh achieved his celebrated result on the basis of a few number of hypothesis, thereby 
capturing the essence of the physics contained in the collapsing cavity. First, he considered 
the fluid to be incompressible, that is, to have a divergenceless flow. By adopting a spherical 
coordinate system, concentric with the cavity, this means that the negative radial component 
of the velocity field v (r, t) may be expressed as 

R 2 R 



(1) 



where R (t) denotes the instantaneous radius of the cavity, and the dot represents total time 
derivative. Further, by assuming a uniform fluid density p = = constant, and applying 
the work-kinetic energy theorem, he was able to relate the collapsing rapidity of the cavity's 
boundary to its instantaneous radius through 



3poo \R 3 

where a is the initial radius of the cavity and stands for the constant fluid pressure 
at infinity. In the above expression, the instantaneous pressure on the cavity's boundary is 
presumed to vanish, once effects provoked by surface tension and viscous terms are neglected. 
Finally, in order to determine the total collapsing time tr of the cavity, Rayleigh integrated 
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Eq. d2J) with respect to R from the initial radius to zero, thereby obtaining 



tr = a 



37TPQQ T (5/6) 

2 PoD r(i/3) 



(3) 



where T (m/n) is the gamma function. Here, we notice that Eq. (j2J) is easily shown to be a 
particular first integral of 



which, from now on, will be referred to as Rayleigh's equation of motion. 

In what follows, we do not try to include any additional phenomenological effect into 
Rayleigh's bubble dynamics, as the ones mentioned in the above quoted references. Instead, 
we seek fundamental modifications of his results by relaxing the hypothesis of uniform density 
assumed by Rayleigh, as well as in some of their important developments like in the so-called 
Rayleigh-Plesset equation [jjj]. 

The article is organized as follows. In section II, we derive an analytical, non-trivial solu- 
tion to the continuity equation, which is consistent with the hypothesis of a divergenceless 
flow. In section III, the A-factor is defined, leading to an extension of Rayleigh's equation 
of motion and a correction to the total collapsing time of the cavity. In sections IV and V, 
we discuss the consistency of the corrections with the Bernoulli theorem and some physi- 
cal consequences in the framework of the so-called Rayleigh-Plesset equation, respectively. 
Finally, in the conclusion section, the main results are summarized. 

II. MASS-SHELL SOLUTION 

For a spherically symmetric one-dimensional non-steady flow, the fluid density p (r, t) is 
demanded to satisfy the continuity equation in the form 



Now, once Eq. ([I]) is assumed to hold true, and changing coordinates from time t to the 
instantaneous radius R(t), Eq. ([5]) readily leads to 



rr + Ir 2 + ^ = o 



2 Poo 



(4) 
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(6) 
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As it appears, one is tempted to take p = constant as the solution to the above equation which 
is consistent with the divergenceless flow. However, by considering the region a < r < b of 
the fluid, a more general solution to that equation may be written as 



, />> Pb Pa e Pa 

P ( r » R ) = — + T—Zx ex P 



3 d3~ 

'' ' (7) 



a 3 a 3 



1 - e" A 1 - e" 

where we have chosen p a = p (a, a) and pi, = p (b, a) as boundary conditions at the initial 
instant, whilst the dimensionless boundary factor A has been defined as 

b 3 

A = ? -l. (8) 

Note that Eq. ([7]) actually describes a vast collection of radially stratified mass-shells at 
some fixed instant. Notwithstanding, since we are describing a finite class of solutions, 
Rayleigh's one may be contained into this as a particular case. Indeed, by identifying a to 
the initial radius of the cavity, and taking the limit b — > oo, from Eq. flTj), it directly follows 
that 

'r 3 R 3S ~ 



p (r, R) = poo + (p - Poo) exp 



a 3 a 3 



(9) 



where we have renamed p = p a and = Pb-^oo- In that case, clearly, p denotes the instan- 
taneous density on the cavity boundary and p^ represents the constant density at infinity. 
Of course, by putting p = poo = constant, Eq. ([9]) trivially recovers the solution adopted by 
Rayleigh. In principle, when a spherical bubble is formed in an uniform medium, one may 
expect a variation of the density with the radial coordinate. Since the mass is conserved 
and a hole has been somewhat created, the mass must be redistributed (compressed) close 
to the frontier of the cavity. In this way, the density near the cavity boundary becomes 
greater than far from the wall, and, solution (j^J), says that such a decay is actually very fast 
(oc e for a divergenceless flow. 

Here, for further purposes, we rewrite Eq. (Q in the form 

p(x) = Poo (l + Ae~ x ), (10) 

where we have changed coordinates from radius r to the dimensionless quantity x through 
the transformation 



for some fixed R (t), as well as defined the dimensionless A-factor by the relation 

A = -^--l. (12) 

Poo 

Given the above analytical possibility of a non-uniform fluid density, as the only new 
modification with respect to Rayleigh's original formulation, and since the assumption of a 
vanishing instantaneous pressure on the cavity's boundary still holds, one can anticipate that 
any correction to the total collapsing time of the cavity must involve only some dimensionless 
combination of both po and p^. Let us now carefully examine the proposed generalization of 
the original Rayleigh equation of motion, which shall show consistency with our non-uniform 
density solution, as well as lead to a correction to the standard total collapsing time of the 
cavity. 



III. GENERALIZED RAYLEIGH'S EQUATION: THE A-FACTOR 

The fluid pressure p (r, t) is required to satisfy Euler's equation in the form 

dv dv 1 dp q 
dt dr p dr 

for a spherically symmetric one-dimensional non-steady flow. Once Eqs. ([T]) and (jlOp are 
assumed to hold true and taking into account Eq. (fTTl) . it proves convenient to introduce 
the dimensionless functional 

y(x) = ^, (14) 

Poo 

for some fixed instant t, or, equivalently, for some fixed instantaneous radius R (t) . Now, by 
integrating Eq. (1131) from the instantaneous boundary of the cavity to infinity, provided the 
instantaneous pressure on the cavity's boundary vanishes, we obtain the extended Rayleigh 
equation of motion 

RR + |# + = 0, (15) 

2 Pc 



oo 



where, from Eqs. ffTUl) and f[T4|) . the time independent dimensionless A-factor is defined by 

A= f_J_§^. (16) 

J 1 + Ae~ x dx v ; 
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As one may check, by applying the initial condition on the surface of the empty cavity 
[R (0) = for R = a], a first integral of Eq. f|T5|) reads 



2A Poo ( 



3 



*=^\*- l Y (17) 

As a result, the corrected total collapsing time of the cavity can be written as 

r A = (18) 

which could be anticipated from the equation of motion (e.g., by replacing p^ — > Ap^). 

All the above results should be compared with Eqs. (jSJ), ([3]) and (jl]) of the introduction. 
In the limit A — > 0, that is, by putting p = p^ = constant, we see that Eq. (fTBj) implies 
A = 1, thereby recovering the respective Rayleigh results. Naturally, in order to quantify 
an arbitrary A-contribution, the integral given by Eq. ( JT6l) must somewhat be evaluated. 
Notwithstanding, we discuss now a simple estimate of A(p /poo), by assuming that A is 
sufficiently large, that is, p Poo- An independent derivation based on the Bernoulli 
equation can be seen at the end of the next section. 

To begin with, we observe that the functionals p (x) and p (x) over the domain of inte- 
gration (0 < x < oo) suggest the definition of a large dimensionless quantity x*, defined as 
the value x takes such that 

y{x*) = \- (is) 

Next, we separate the integral in Eq. f|T6l) in two parts, namely, the first over the interval 
< x < x*, and the second one for < x < oo. Since A is assumed to be large enough and 
the cavity very small as compared to the infinite medium, we may neglect the functional 
dependence of the integrating factor on x throughout the first interval, whilst keep it through 
the second one. However, since x is also large throughout the second interval, the whole 
integrating factor may be fairly well approximated to unity over there. As a consequence, 
we may approximate our Eq. (jT6l) by 



from which it immediately follows that 

= EL 

2p 



A = P0 + P °°. (21) 
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Note that, as long as the instantaneous density on the cavity's boundary severely exceeds 
the constant density at infinity, A approaches its lowest possible value, namely 

A m i n = — . (22) 

As a consequence, the corrected total collapsing time of the cavity approaches its largest 
possible value, namely 

W = t r V2. (23) 

Let us now investigate the consistency of the A-factor as given by Eqs. ffl6|) and (|2ip 
with Bernoulli's theorem. As we shall see, the foregoing independent argument leads to the 
same A-factor and Rayleigh's type equation of motion. 



IV. A-FACTOR AND THE BERNOULLI THEOREM 

The fluid motion represented by Eq. (pQ) implies V X v — 0, i. e., the velocity field in 
the Rayleigh dynamics yields a potential flow. As a consequence, we can immediately write 
v = V0, where the hydrodynamical potential <j) (r, t) may be chosen to be 

R 2 R , , . 

</> = . 24 

r 

Now, for an isentropic flow, viz., on the assumption of a uniform specific (per unit mass) 
entropy s throughout the fluid, a generalized form of Bernoulli's theorem shall be written 



as 



m 



f + \v 2 + u = f{t) 1 (25) 

where u (r, t) denotes the fluid enthalpy, and / (t) represents some arbitrary function of 
time. At this point, we can, simultaneously, evaluate the left-hand side of Eq. (|2"5"|) on the 
cavity's boundary, that is, at r = R, and, consistently with Eqs. ([T]) and (1241) . which ensure 
that both v (r, t) and (r, t) vanish in the limit r — > oo, choose / (t) to be the constant value 
lu (r, t) takes at infinity. In that case, once uj (r, t) is related to the specific internal energy 
e (r, t) of the fluid through 

oo = e + ^, (26) 



by assuming a vanishing instantaneous pressure on the cavity's boundary, Eqs. (I25I) and 
(|26|) lead to 



RR + \k 2 + + ( eoo - 6 ) = 0, (27) 

^ Poo 



where eo and e m denote the constant values e (r, t) takes at the cavity's boundary and at 
infinity, respectively. It thus follows that any correction to the original Rayleigh's equation 
of motion may only come from the parenthesis in Eq. ( 1271) . In order to evaluate it, we shall 
make use of the well-known thermodynamical relation 

de = Tds + p^, (28) 
P 2 

where T (r, t) denotes the thermodynamical temperature throughout the fluid. However, 
since we have an isentropic flow (ds = 0), we can rewrite Eq. (|27|) as 

RR + \k 2 + — + r V % = 0. (29) 
2 Poo J po P 2 

It is a trivial matter to show that the above equation is equivalent to the previously derived 
form of the generalized Rayleigh's equation of motion, Eq. ffl5l) . provided the integral in 
Eq. (1291) is carried out by parts, 

px d p poo r id p, 

p— = + / -Trdr. (30) 

po P Poo Jr par 

As a consequence, from Eq. (1271) . we arrive at an independent definition of the A-factor 
with respect to the one given in Eq. ffl6l) . 

A = l + ^( eoo - eo ) = l + ^ V%- (31) 

Poo Poo J P0 P 2 

At this point, through a quite simple argument, consistently with a non- uniform den- 
sity solution, we show that the same estimate of the A-factor as given by Eq. (|2T|) can be 
achieved in the framework of Bernoulli's theorem. The unique necessary assumption is the 
monotonical growth of the fluid pressure throughout the increasing radial coordinate con- 
centric with the cavity. We start by introducing an instantaneously averaged fluid pressure, 
p say, defined throughout the bulk of the infinite medium as 

p= lim ^ [Sp + 25p + 35p + . . . + (N — 1) Sp + p^] , (32) 



s 



where N denotes a very large number of mass-shells, on the assumption that p (p) on each 
radially stratified layer. In principle, one should expect p to be dependent on the fluid en- 
tropy, as well. However, we have assumed the flow to be isentropic. In this proposed scheme, 
two neighboring surfaces, no matter if they stand close together or far apart throughout the 
radial direction, are labeled by two different values of p (p), whose difference is exactly equal 
to a very small number, Sp, say. Note that p^ = N5p, whilst the pressure on the cavity's 
boundary presumably vanishes. In addition, consistently with the assumption that p^ = 
constant, 5p — > in the limit N — > oo. Now, the problem of evaluating the sum in Eq. (|32|) 
is a trivial one. The answer is 

5p + 25p + 35p+... + (N-l)5p + Poo = N ( 6 P + P°°\ (33) 

In other words, the sought instantaneously averaged fluid pressure, actually, does not depend 
on time, it is a constant, and we can legitimately approximate the pressure throughout the 
volume of the fluid by 

P= P ~f- (34) 

As one may check, by inserting Eq. into Eq. ( I3T1) . our previous expression for A as 
given by Eq. (1211) is recovered. 

V. A-FACTOR AND RAYLEIGH-PLESSET EQUATION 

It is also interesting to discuss how the A-factor (or the variation of the density) may 
modify the description of other fluid properties, from the point of view of bubble dynamics 
theory. In this regard, an important analytical development was achieved by Plesset thus 
leading to what is now often referred to as Rayleigh-Plesset equation. That extended theory 
is remarkably simple and, in particular, explains several features related to single-bubble 
sonoluminescence experiments. In such an approach, the fluid density is also constant but 
one assumes a gas-filled cavity so that several boundary effects are taken into account as 
well. The modified equation of motion reads 1(3, 12] 



RR+1r* + p °° p(jR) =0, (35) 

2 poo 
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where p(R) is the pressure of the liquid at the bubble boundary. It includes terms like 
surface-tension, viscous stress and the pressure in the bubble, pb, and can be written as 

p(R)=p B ---4^, (36) 

where a is the surface-tension constant and p is the viscosity coefficient. From the above 
expression, the Rayleigh-Plesset equation assumes the following form 

RR+ |# + J_ L + ^ + 4 p| ) = 0. (37) 
Z Poo \ K K I 



An important observable difference between Eqs. (1371) and (J4j) is the existence of an equili- 
brium (unstable) radius, -Re, which is obtained by imposing the conditions R = and R = 
to be simultaneously fulfilled. For further reference, it is given by 

Re = — ^— , (38) 
Pb - Poo 

which does not depend explicitly on the fluid density. Additional properties of the Rayleigh- 
Plesset equation have been widely discussed in the literature [id. 111]. 

In the above context, the interesting question is: how the non-uniform density solution 
of section II will affect the Rayleigh-Plesset equation? Naturally, one may also expect a 
modification of the standard equilibrium radius. 

Following the same approach presented in Section III, it is readily checked that the 
Rayleigh-Plesset relation, Eq. (I3"7|) . must be replaced by 

HB + *# + L*>> =0 . (39) 
2 poo p R po R po 

As one should expect, beyond the modification of the A-factor, all boundary terms are now 
divided by the density of the liquid in the evolving bubble frontier (p ). In addition, the 
equilibrium conditions lead to a corrected equilibrium radius of the bubble 

(A) = 2a , (40) 

Pb - A (po/Poo)Poo 

which depends explicitly on the boundary value of the density. It is also worth notice that 
in the limit A = 1, that is, for po = Poo = constant, Eq. ( f40l) trivially recovers the standard 
equilibrium radius of the bubble as predicted by the standard Rayleigh-Plesset equation. 
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VI. CONCLUSION 



In this work, we have discussed a new kind of contribution to Rayleigh's equation of 
motion and the total collapsing time of an empty spherical cavity. As we have seen, a 
correction term, herein called the A-factor, naturally appears if the fluid density is not 
uniform, as assumed by many authors including Rayleigh himself. Indeed, the A-factor is 
a direct consequence of a divergenceless spherically symmetric non-steady flow, whose fluid 
density is an evolving function in space and time, as given by the general solution derived 
in section II [see Eq. 

The inclusion of the A-factor as a new parameter into Rayleigh's dynamics necessarily 
leads to non-negligible corrections of a number of meaningful physical quantities, as the 
collapsing time. We have also shown how the extension usually named Rayleigh-Plesset 
equation is affected by the non-uniform time varying density. In particular, a new equi- 
librium radius has also been derived. Unlike the standard approach, it depends explicitly 
on the fluid density at the boundary and far from the bubble. All the results summarized 
above are analytical and have been justified from first principles. It should also be stressed 
that even more general formulations, as the one describing approximately the influence of 
the liquid compressibility 15], may be affected by the correction proposed here. A more 
detailed investigation of the A-factor on bubble dynamics will be discussed in a forthcoming 
communication. 
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